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Abstract. This paper presents observations of an extended Lyman-a emission nebula surrounding the galaxy responsible for 
the Damped Lyman-a Absorption (DLA) line in the spectrum of the quasar Q2233+131. With the Potsdam Multi Aperture 
Spectrophotometer (PMAS) we measure the properties of the extended Lya emission in an area of 3" x 5" having a total line 
flux of (2.8 ± 0.3) X 10^^*^ erg cm"^ s"\ which at redshift z = 3.15 corresponds to a luminosity of (2.4t° J) x 10*^ erg s"^ 
and a size of 23 x 38 kpc. The location of the emission is spatially coincident with the previously detected DLA galaxy, but 
extends significantly beyond its limb. We argue that the Lya emission is likely to be caused by an outflow from the DLA 
galaxy, presumably powered by star formation. In the case of negligible dust extinction, the Lya luminosity indicates a star- 
formation rate of 19 ± 10 Mq yr~^ consistent with that derived from the UV continuum flux from the parent galaxy. The wind 
velocity indicated by the integral field spectra is of the order of several hundred km s~^. We find no indication of emission 
originating in a rotating disk. 
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1. Introduction 
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High redshift quasars (QSOs) show multiple absorption lines 
bluewards of the redshifted 1216 A Lymana (Lya) wave- 
length. Known as the Lymana forest, this is caused by ab- 
sorption in neutral hydrogen clouds along the Une of sight to- 
wards the QSO. Clouds having column densities larger than 
2x lO'^^cm^'^ give rise to line profiles with broad wings char- 
acteristic of damped Lya lines. To date approximately 150 
Damped Lya Absorbers (DLAs) with redshifts of 0.1 < 
z < 4.6 have been confirmed dCurran et alJ l2003) . It has 
been found that DLAs contain a significant fraction of to- 
tal gas mass co mpared to the mass of the stars in present 
day galaxies (iWolfe et alJll995HStorrie-Lombar di et al. 1996; 
IStorrie-Lombardi & Wolfel2000l) . and it is questioned whether 
there is a significant evolution with redshift ( Rao & Turnsheld 
[2OOO) . Spectroscopic observations have shown that the DLAs 
have metalUcities of 0.01 -1 times solar with a m ild increase 
with decreasing redshift ('Proc haska et alJ l2003l) . suggesting 
that DLAs are star-forming objects. Nevertheless, the rela- 
tion between DLAs and galaxies is not well understood. At 
higher redshifts, the DLA galaxies have been suggested to be 
thic k disks or the pro genitors of present day spirals galax- 
ies llWolfe et alJll98^ . while others suggest that the counter- 
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parts could be dwarfs (iHunstead et'aDll990l) . o r galaxy build- 
ing blocks in a hierarchical merging scenario (iHaehnelt et alJ 
.1998 ). 

Many investigations have been performed in order to 
establish what the galaxy couiiterpart s to the DLAs re- 
semble most (e.g. 'Le Brun etaP Il997t IWarren et all l200lt 
Colbert & Malkan 2002). Typically, deep broad-band or 
narrow-band imaging of the fields containing the DLAs has 
been carried out, and objects near the line of sight of the QSOs 
are detected after subtraction of the QSO point spread function. 
These candidate DLA galaxies are typically faint. Successive 
follow-up spectroscopy of the candidates is required to reveal 
whether or not they have the same redshift as the DLA line. 
Only in 4 cases DLA galaxies have been c onfirmed this way 
for the high redshift (z >, 1.9) DLA galaxies ('M0ller & W arrenj 
11993; Diorgovski et al. 1996; Fvnboetal. 1999; Miller etall 
|2p02). In two additional cases Lya emission lines have been 
detected in the troughs of the DLA lin es in the Q SO spectra 
(Leibundgut & Robertson 1999; Elliso n et alJ2002l) . 

We here present a study of the Q2233-H131 at z = 3 .295, 
which has a DLA line at z = 3.153 dSargen t et al."l989^. The 
metallicity of the DLA is [Fe/H] = -1.4 ( Euetal. 199_7) , and 
the column density of H I is below the classical limit of a DLA 
line, having A'hi = 1 x 10^" cm~^, thus formally characteriz- 
ing this as a Lyman-Limit system. In accordance with previous 
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papers on this object we will continue to denote it a DLA ab- 
sorber. A candidate galaxy responsible for the absorption was 
found at an impact parameter of 2'.'3 using the Lyman break 
techn ique, suggesting a redshift larger than 3 dSteidel et alJ 
Il995h . This object was confirmed as th e absorbing ga laxy hav- 
ing the same redshift as the DLA line jDiorgovski et alJI 19961 
hereafter D96). These authors measured a Lya line flux of 
(6.4 ± 1.2) X 10^^^ erg cm^^ s^^ in a Keck long-slit spec- 
trum and found the magnitudes of the Ly man break ga l axy of 
R = 24.8 ± 0.1 and V = 25.1 ± 0.2.IWarren et alJ ll200lh 
found H = 25.34 ± 0.17 for the ga laxy usin g near- IR p hotom- 
etry with the HST/NICMOS, and lM0Uer et alJ <2002i) found 
^50 = 25.75 ± 0.12 with STIS images. 

With integral field spectroscopy the conventional two step 
approach for confirming a DLA galaxy can be avoided. A 
previous attempt to use this technique only yielded an upper 
limit for the line emission from the DLA galaxy in BR 1202- 
0725 dPetitiean et al .119961) . while the spectral range of the ob- 
servations of the DLA system in APM 08279H-5255 did not 
cove r the appropriate w avelengths for the redshifted Lya emis- 
sion jLedoux et al .'^1998). We show here that with the Potsdam 
Multi Aperture Spectrophotometer (PMAS) instrument we not 
only detect Lya emission from the DLA galaxy at z = 3.15, 
but we find that the object causing the Lya emission is ex- 
tended, and the line flux is larger than reported previously in 
the literature. 

In Sect. |2] of this paper we will describe the spectroscopic 
observations, and the procedures for reducing the data. We 
analyse the spectra of the QSO and the extended emission from 
the DLA in Sects. |3]and0] respectively, addressing the nature 
of the extended emission. The relation between the Lya emis- 
sion and the location of the DLA galaxy seen in high spatial 
resolution Hubble Space Telescope (HST) images is described 
in Sect.|5] Similarities with other Lya emitting objects are de- 
scribed in Sect. |6] which leads to possible interpretations for 
the origin of the extended Lya nebula in Sect. The ionised 
gas mass in the nebula is estimated in Sect. |8l In Sect. |9] we 
present our conclusions. 

Throughout the paper we assume a flat Universe with Hq = 
70 km s-i Mpc-\ 0,^ = 0.3, and = 0.7. The redshift z = 
3.15 then corresponds to a luminosity distance of 8.3 x 10^^ cm 
(27.0 Gpc), 1" corresponds to a linear size of 7.6 kpc, and the 
look-back time is 11.5 Gyr. 

2. Observations and data reduction 

The PMAS integral field instrument uses two cameras: A cryo- 
genic acquisition and guiding camera (A&G camera) that can 
be used for i maging in add ition to the integral field spectro- 
graph (IFS) ('Rot h et allEob O). The A&G camera has a SITe 
TK1024 chip with Ikx Ik pixels with a scale of 0'.'2 per pixel 
giving a field of view of 3f4x3.'4. The PMAS spectrograph is 
equipped with 256 fibers coupled to a 16x16 lens array, that 
we used with a spatial sampling of 0'.'5xO'.'5 per fiber on the 
sky, resulting in a field of view of 8"x8". The spectrograph 
camera has a SITe ST002A 2kx4k CCD and its 256 spectra 
have a FWHM of ~2 pixels when using a 2 x 2 binned read-out 



mode. The spectra are aligned on the CCD with 7 pixels be- 
tween adjacent spectra making cross-contamination negligible. 

PMAS is mounted on the 3.5m telescope at Calar Alto. 
We observed Q2233H-131 on Sep. 2 2002 for a total of 7200 s 
(4x 1800 s) at an airmass between 1.094 and 1.15 and a seeing 
between I'.'O and 1'.'3 measured by the A&G camera. Using a 
300 gr/mm grating yielded a spectral resolution of 6 A while 
the grating was set to cover the wavelength range of 3930- 
7250 A. Calibration images were obtained following the sci- 
ence exposures and consisted of spectra of emission line lamps 
(HgNe), and spectra of a continuum lamp needed to locate the 
256 individual spectra on the CCD. Observations of the spec- 
trophotometric standard stars BD H-28°421 1 and Hz4 were ob- 
tained during the night for flux calibration. 

Reduction of the data was done in IDL with P3D_online, a 
softw are pack age written specifically for reducing PMAS data 
(Bec keJ2002h . After bias subtraction the 256 spectra were ex- 
tracted from the two-dimensional frames. Wavelength calibra- 
tion was performed using the emission line lamp spectra. For 
the flat fielding the average transmission of each fiber was de- 
termined using exposures of the sky at twilight. Cosmic ray hits 
were removed from each of the 4 files u sing the L. A. Cosmic 
routine within IRAF ( Ivan Dokkuml200ll) . and the rejected pix- 
els were inspected by eye checking that no pixels close to the 
Lya emission line were affected. The data cubes were cor- 
rected for the effe ct of differential a tmospheric refraction us- 
ing the formula of Filippenko ( 1982). Given the small airmass 
of Q2233-I-131 during the observation the effect is small, but 
not negligible at long wavelength intervals. At the wavelength 
of the DLA line the effect can be ignored, since this was the 
reference wavelength for estimating the relative offsets before 
combining the 4 frames. This method ensures that the posi- 
tion of the DLA galaxy relative to the centroid of the QSO at 
other wavelengths is not shifted. In the end the four files were 
co-added, resulting in a data cube of dimensions 16 x 16 x 1024 
pixels. 

For subtraction of the sky background, an average back- 
ground spectrum was created by co-adding several spectra at 
the edge of the field of view, uncontaminated by the QSO 
flux, and subtracted from all 256 individual spectra. All further 
data calibration was done using IRAF. For the standard star 
observations we coadded all spectra within a radial aperture 
of 3", and compared the one-dimensional standard star spec- 
trum with table values to create a sensitivity function taking 
into account the atmospheric extinction typical for Calar Alto 
( faopp & Fernandez. 2QQ2) . Finally, the spectra were flux cali- 
brated using this sensitivity function. 

Further analysis and inspection of the data was performed 
with the "Euro3D Visualization Tool", whic h is a very effic ient 
tool made for visuaUzing integral field data (ISanchezl2003h . 

3. Tlie QSO spectrum 

We constructed narrow band images by selecting appropriate 
wavelength intervals in the combined 3D data cube, resulting in 
images with 16x16 spatial pixels ("spaxels"). In Fig.^such an 
image of Q2233-I-131 is shown in the wavelength range 5200- 
5300 A. For creating a one-dimensional object spectrum one 
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Fig. 1. Left hand panel: Image of Q2233+131 in the wavelength range 5200-5300 A. This image shows the 16x16 spatial pixel 
("spaxel") field of view of PMAS corresponding to 8" x 8". The QSO is centered in the field at (0,0). Right hand panel: Co added 
narrow band image of the Q2233+131, with a field of view of 8"x8", north is up and east is left. Selecting the wavelength 
range 5040-5055 A an extended object of roughly 5"x3" appears, which is the extended Lya emission from the DLA galaxy. 
The selected wavelengths correspond to the DLA absorption trough in the QSO spectrum. The centroid of the QSO is at (0,0) 
while the position of the DLA galaxy found from the broad band observations in D96 and in the space-based images is at (- 
0'.'9,-2'.'2). Apparently, the extended source overlaps with the position of the QSO, but this is caused by the fact that some QSO 
emission is present at the longest selected wavelengths (A > 5050A), as emission from the red wing of the DLA line has been 
included. Contours of 5-8(7 levels above the background have been overplotted for guiding the eye. The contours correspond to 
the combined emission from the Lya nebula and the QSO. The high signal at the left edge of the field at (3.5,1) is caused by a 
bad fiber. 



can select spaxels in this image, each of which represents a 
single spectrum, and co-add the selected spectra. 

For comparison with later spectra we created a combined 
spectrum of Q2233H-131 by co-adding all spectra within a ra- 
dial aperture of 2" (~45 spaxels) in Fig.^ An overall aperture 
correction of 2% was applied, and the resulting spectrum is 
shown in Fig.|2l 

We detect many absorption features from metal lines in 
the spectrum of the QSO associated with the z — 3.15 DLA 
system, some of which were already recognized in D96. The 
redshifts of these lines, listed in Table Q matches the red- 
shift of the DLA line. We derive a mean systemic redshift 
z = 3.1475 ± 0.0005 from the all metal lines in the DLA apart 
from the Si IV and O IV 1400 A blend which has a lai-ger off- 
set than the low ionization species. A detailed analysis of the 
metallicity and column densities is outside the scope of this 
paper, as higher spectral resolution would be necessary for this 
purpose. For the same reason the Ly/3 absorption line is blended 
with the Lymana forest. The dip seen in the tip on the Lya 
emission from the QSO is probably caused by Si II 1260 A at 
the redshift of the DLA system. 

In addition to the absorption lines we find broad emission 
lines from Lya, N V 1240, Si IV H- O IV 1400 A and C IV 1549 
A from the QSO. Including these four lines we find the redshift 
z = 3.2877 ± 0.0052 for the QSO. 



Ion 


Alab [A] 






Lya 


1215.67 


3.2949 




NY 


1238.82, 1242.80 


3.2840 




Si IV + IV] (blends) 


1400 


3.2881 




CIV 


1548.20, 1550.78 


3.2838 




Lya 


1215.67 




3.1476 


Si II 


1260.42 




3.1468 


OI 


1302.16 




3.1479 


CII 


1334.53 




3.1476 


Si IV 


1400 




3.1490 


Sill 


1526.71 




3.1485 


CIV 


1548.20, 1550.78 




3.1477 


AlII 


1670.789 




3.1478 


Lya(DLA galaxy) 


1215.67 


3.1538 





Table 1. List of emission lines from the Q2233H-131 itself and 
the absorption lines, related to the z — 3.15 DLA system 
detected in the QSO spectrum. A standard air-to-vacuum cor- 
rection has been applied to the observed lines before deriving 
the redshifts. The uncertainties of the measured redshifts are 
typically ±0.0005. 



4. The spectrum of the DLA galaxy 

In Fig. [Owe show another narrow-band image, this time se- 
lected from the wavelength interval 5040-5055 A, correspond- 
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4000 4500 5000 5500 ^ 6000 6500 7000 

Wavelength [A] 

Fig. 2. Spectrum of Q2233+131. All spaxels within a radial aperture of 2" have been coadded and then scaled to the total flux 
of the QSO, measured within a 3" radial aperture. The bright line at 5200 A is Lya emission from the QSO. This and other 
emission lines from the QSO are labeled under the spectrum. The DLA line at 5050 A and the corresponding metal absorption 
Unes belonging to this system are labeled above the spectrum. Redshifts of all lines are given in Tabled 



ing to the absorption trough in the QSO spectrum. An extended 
object is visible to the south and south-west of the location of 
the QSO. The size of this object is roughly 5"x3", however, 
some of the extended object emission is due to the presence of 
emission in the red wing of the DLA line from the QSO. We 
note that the very high signal which appears at the left edge of 
the field with coordinates (3.5,1) is due to a bad flat-field effect 
of one single fiber, which has a lower overall transmission. 

We have added 35 spaxels which are apparently associated 
with the emission line object. A part of the spectrum around 
the emission line is shown in the upper panel in Fig. |3] where 
one clearly sees an emission line at ~5050 A. This wave- 
length corresponds to Lya at the redshift of the DLA line at 
z — 3.1538 ± 0.0005, in agreement within la with the redshift 
published in D96. A second line is present at ^5200 A which 
is caused by the Lya emission from the QSO at a distance of 
2". For comparison the same section of the QSO spectrum is 
shown in the middle panel in Fig.|3] There appears to be another 
broad absorption line at 5000 A, but a high resolution spectrum 
of the QSO has shown that this f eature is caused by a blend 
of 5-6 individual absorption lines dSechtoldl 19941) . The lower 
panel shows the spectrum of the DLA where a scaled spectrum 
of the QSO has been subtracted. In fact, the Lya emission line 
at 5050 A is effected very little by this subtraction. 

The spectra were analysed with the ONEDSPEC package 
in IRAF. We measured a full width half maximum (FWHM) 
of the Lya emission fine of 20 ± 2A. The resolution in the 



combined spectra is 7.9 A, measured from the FWHM of the 
5577 A night sky line in the combined spectra. This yields an 
internal FWHM of 18.4 ± 3.3 A of the Lya fine, coiTespond- 
ing to a restframe velocity of 1090 ± 190 km s^^. We see no 
emission Une from N V A1240 or C IV A 1549 from the DLA 
galaxy down to a 3cr detection limit of 1 x 10^^^ erg cm^'^ s^^ 
in agreement with D96, albeit their detection limit was 10 times 
fainter. 

The redshift difference between the DLA galaxy's Lya 
emission component and the DLA absorption components in 
the QSO spectrum is Az = 0.0063±0.0007 corresponding to a 
velocity of 450±50 km s^^. This is larger than the 209 km s^^ 
reported in D96. Here we have used information from all the 
detected absorption lines in calculating the systemic redshift of 
the DLA. If we instead used the redshift reported in D96 for 
the absorption minimum of the DLA line itself, we would find 
a velocity difference of 270±40 km s^^, which is within 2a 
of their value. Since the interstellar absorption lines have been 
shown to be blue-shifted by > —100 km s~^ with respect to the 
stars in Lyman Break Galaxies jShaplev et al.-2003.) . the latter 
value (^300 km s^^) is probably a better estimate. 

The line flux measured by fitting the observed line with 
a Gaussian profile is (2.4 ± 0.3) x 10^^^ erg cm^^ s^^; a 
simple summation of the flux values for each pixel from the 
emission line give s the same result wit hin the errors. From 
the dust maps of ISchlesel et al.' ('l998') a Galactic redden- 
ing of Eb-v = 0.068 in the direction towards Q2233H-131 
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Fig. 3. Spectrum of the extended emission associated with the 
DLA galaxy at z = 3.15 in Q2233+131. The spectrum has 
been obtained by co-adding ~35 spectra which appear to be as- 
sociated with the DLA galaxy in the image in Fig.^ The few 
spaxels which are overlapping the position of the QSO have 
not been coadded, since these are contaminated by flux from 
the wing of the QSO and not associated with the Lya emission 
object. An emission line at 5050 A corresponds to Lya at the 
redshift of the DLA line in the QSO. The broad emission fea- 
ture at ~5200 A is due to residual Lya emission from the QSO 
at a radial distance of 2" from the centre of the QSO. For com- 
parison we have plotted the total QSO spectrum in the middle 
panel with the same wavelength range. Note the scale on the 
y-axis differs by a factor of ~30 between the two plots. In the 
lower panel we have plotted the DLA spectrum where a scaled 
QSO specti'um has been subtracted. 



is found. Correcting for this effect increases the line flux to 
(2.8 ± 0.3) X 10^^^ erg cm^^ s^^ and the total luminosity of 
the source is 2.4^q'2 x 10^"^ erg s^^ in the adopted cosmol- 
ogy. With the noise in each spectral element around 5000 A 
of 3.2 X 10^^* erg cm^-^ s^^ A^^ and the line flux being mea- 
sured over 20 pixels, the significance of the detection of the 
Lya line is 15a. The size of the object is 23x38 kpc above a 
3a detection threshold of 1 x 10^^^ erg cm^'^ s^^ A^^. 

A measurement of the equivalent width (EW) of the emis- 
sion line is hampered by the fact that the continuum emis- 
sion is extremely faint. From the observed ground-based V 
and i? band magnitudes in D96 together with magnitudes from 
HST images (derived in Section O, we estimate the flux from 
the underlying continuum at ~5050 A to be (LI ± 0.4) x 
10^^^ erg cm^^ s^^ A^^ assuming a power-law spectral en- 
ergy distribution, fi, oc i^^ in the continuum. With this line and 



continuum flux we derive the rest frame EW=190j^7Q" A. The 
large error is mainly caused by the uncertainty for the contin- 
uum flux. Furthermore, the measured EW must be taken as a 
lower limit because the effect of dust extinction affects the Lya 
flux more than the UV continuum emission. 



4. /. Velocity structure 



With integral field spectra one can in principle determine the 
velocity structure of the Lya nebula. Specifically, it would be 
interesting to search for signs of rotation. A detailed analy- 
sis of the velocity structure is not possible with the present 
dataset, since in several of the individual spectra the emission 
line is only detected on a r^2a level, which makes any anal- 
ysis of any velocity structure unreliable. Instead we summed 
up the ~15 spectra in the left-most (eastern) part and the 
^15 in the right-most (western) part, respectively. A sketch 
of the division is shown in Fig. |3 The Lya emission for 
each region is still detectable, and the line flux of each is 
~lxlO^^^ erg cm^^ s^^ A^^. We fitted Gaussians to the 
emission line in these two regions, and find that the peak of 
the emission is shifted by H-2.5 A from the east to the west 
part, which corresponds to a rest-frame velocity difference of 
^^150 km s^^. This is small compared to the rotational veloc- 
ity of present day large spiral galaxies. Repeating this exercise 
by splitting the extended object into a southern and a north- 
ern part as shown in Fig. 0] the emission line is shifted by 5 
A corresponding to a difference of ^300 km s~^, where the 
northern part has the largest redshift. This velocity is compa- 
rable to that expected for spiral galaxies, but the geometry is 
unusual for a disk, and as shown in Sect. I4.2l the apparent shift 
of the emission line can be caused by a combination of more 
than one emission region. In Table |2] we summarize the prop- 
erties for the different parts of the extended Lya nebula, and in 
Fig. Is] we show the Lya line for the different regions. The la 
errors for the centroids of the Gaussian profiles are 0.5 A esti- 
mated from simulations of artificial spectra having faint emis- 
sion lines at known wavelengths. From these simulated spectra 
we also estimate that the la error for the EWs is ~2 A. All 
regions are seen to have a large FWHM indicating velocities of 
700-1 100 km s^^, i.e. they are well resolved by our spectra. 

Previously, the existence of a disk in this system has 
been suggested from studies of the metal absorption line pro- 
files ( Lu et al. 1997 ). This interpretation is not supported by 
iLedoux et alJ ( fl998li) . who concluded that the profiles could 
be caused by several interacting components, an interpreta- 
tion which is motivated by cold dark m atter simulations of hi - 
erarchical clustering of proto-galaxies dHaehnelt et aljfl998i) . 
Under the assumption of a rotating disk, the enclosed mass cal- 
culated from A/dyn = v^r/G is ^ 10^'' Mq. As will be dis- 
cussed below, the Lya emission is unlikely to originate in a 
rotating disk, and the mass calculated here is probably overes- 
timated. 
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Fig. 4. The images show the lowest (5a) contour level taken from Fig. [2 The chosen division between a north and a south region 
is indicated in the left panel, and the division between an east and a west region to the right. In both cases the emission spatially 
coincident with the QSO has not been included. 



region 


centroid 


/line 




FWHM 


EW 




(A) 


( erg cm 


"^ s-^) 


(kms"^) 


(A) 


total 


5048.3 


2.4x10" 


16 


1090±190 




east 


5049.6 


1.0x10" 


16 


770±150 


west 


5052.1 


1.3x10" 


16 


1120±200 




south 


5048.8 


0.7x10" 


16 


650±150 




north 


5053.4 


1.3x10" 


16 


890±150 





Table 2. Properties of the Lya emission from different parts of 
the extended nebula derived from fitting a single Gaussian pro- 
file to the observed spectrum. The la errors of the estimations 
of the Gaussian centroids are 0.5 A. Note that there is an over- 
lap of the selected spectra belonging to the various regions such 
that some spectra belonging to the east part also belong to the 
north (see Fig. 1^. The EW is not meaningful to calculate for 
the different regions, as the continuum emission is restricted to 
the total emission. Line fluxes listed in column 3 have not been 
corrected for Galactic extinction. 




5000 5020 5040 5060 
Wavelength [A] 



5080 



5100 



Fig. 5. Spectra of the four different regions of the emission line 
object. The north, south, east, and west regions correspond to 
those listed in Table |2 All regions are seen to have the same 
double peaked profile of the Lya line. 



4.2. The double peaked Lya line 

In Fig. 12] the Lya emission line appears to have a double 
peaked profile. As this feature is not apparent in D96, we dou- 
ble checked carefully for its integrity. Although the dip near 
the line centre is of the order of the noise level, it appears in 
all spectra and therefore is almost certainly real. In Fig. [S] all 
the four separate spectra described above have the same double 
peaked profile. The feature is not caused by the subtraction of 
the sky background, nor is it due to imperfections in the CCD 
in the region around the emission fine. 



The difference between our measured emission redshift and 
the one reported in D96 may indicate that our wavelength cali- 
bration has an error of ^1 A. This would not have an effect on 
the reported velocity difference since the error would be sys- 
tematic in all our spectra and not affect the differences in the 
estimated Gaussian maxima. Furthermore, as mentioned, the 
spectra are contaminated by the QSO, but a subtraction of a 
scaled QSO spectrum from the four individual spectra does not 
change the Lya line profile. 

Anyhow, we caution the reader that other features appear 
to be systematic in the spectra, e.g. a slightly smaller system- 
atic depression is present in all the spectra around 5020 A that 
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region 


centroid 1 
A 


centroid 1 
A 


abs. centroid 
A 


east 


5045.9 


5056.3 


5051.1 


west 


5044.8 


5057.3 


5051.0 


south 


5046.3 


5056.3 


5050.7 


north 


5046.5 


5058.8 


5052.2 



Table 3. The location of the two peaks from fitting simuha- 
neously two Gaussians to the spectra in Fig. [S] The centroids 
for the fits between the 4 regions vary with < 2.5 A indicating 
small velocity differences. For all regions the splitting between 
the two Gaussians are 10-12 A corresponding to velocities of 
600-750 km s~^. The third column lists the centroid for the 
absorption line for the fits consisting of an emission line and an 
absorption line. 



could be due to an error in the background estimation. Keeping 
this in mind, we note that a similar double peaked feature was 
observed for the Lya line of the DLA galaxy in Q2059-360. 
How ever, for this object the hne profi le changed with slit posi- 
tion jLeibundgut & RobertsonI 19991) . 

The values for the centroids of the Gaussian fits listed in 
Table|2]can thus be affected by unequal contributions from two 
separate emission components. We re-analysed the profiles by 
simultaneously fitting two Gaussians resulting in centroids of 
the two peaks listed in column 1 and 2 in Table |3] For these 
two-Gaussian fits the widths of the two individual components 
are barely resolved with this spectral resolution. The location 
of the first peak is shifted by 1.1 A from east to west, corre- 
sponding to 65km s~^, while the second component is shifted 
by-lA. Since the uncertainties for the estimate of the Gaussian 
centroids are 0.5 A, the shifts are consistent with within la- 
errors. The shifts for the north-south regions are likewise small 
(< 2.5A corresponding to < 150 km s^^). For all regions the 
splitting between the two Gaussians are 10-12.5 A correspond- 
ing to velocities of 600-750 km s~^. Very likely, the velocity 
structure across the extended Lya nebula is complex, and sin- 
gle Gaussian fits for such extended areas as analysed here is 
over-simplistic. 

Additionally, we have fit the profiles by a combination of 
a single Gaussian emission line and a Gaussian absorption line 
creating the dip seen in the spectra. The centroids for the ab- 
sorption line for these fits are listed in Table |3l in column 3. 
The x^s from these fits were slightly smaller than for the two 
emission component fits, indicating that the presence of an ab- 
sorption line is preferred. 

4.3. Artificial slit spectra 

With integral field data we can reproduce the observations 
expected from slit spectroscopy. In Fig |6l we show a cut of 
the PMAS spectra around the DLA line in the QSO, where 
we have indicated the emission from the extended Lya neb- 
ula by small circles. In comparison we have in the upper 
right hand panel created an artificial long sht spectrum with a 



width of 1", and the corresponding one-dimensional spectrum 
is shown in the lower panel. These plots can directly be com- 
pared to the ones presented in D96. We derive a line flux of 
(6.5 ± 1.0) X 10^^'' erg cm^^ s^^ A"^ for the Lya line, con- 
firming the flux reported in D96, and a FWHM of 6.5 ± 1.0 A, 
corrected for the instrumental resolution, corresponding to a 
velocity of 390 ± 60 km s^^. This FWHM is similar to the 
the value in D96. With these simple exercises we have shown 
that it is possible to reproduce the previously published results 
derived from long-slit spectroscopy. 

One sees furthermore in Fig. |6l that the Lya emission in- 
dicated by the small circles do not overlap spatially with the 
location of the QSO. This implies that the contours in Fig.[2at 
a distance less than 1" from the QSO only show emission from 
the QSO itself and not Lya emission from the DLA galaxy. 
Therefore we consider that selecting the emission region from 
the DLA galaxy as indicated in Fig.0]is the real extension of 
the Lya nebula associated with the DLA galaxy. 

5. HST deep imaging 

For comparison of the extended emission with high spa- 
tial resolution and deeper optical data previously unpublished 
HSTAVFPC2 images of Q2233H-131 were retrieved from the 
HST archive. Previously, a high s patial resolution image from 
HST/STIS was presented by Meal ier et all 1I2OO2I) and we also 
retrieved the STIS images from the archive for comparison. 
The WFPC2 planetary camera images were obtained through 
the F702W filter with a total integration time of 10800 s and 
the images were combined using the drizzle package in IRAF 
(iFruchter & Hook..2002) . The resulting pixel scale in the driz- 
zled image is 0'.'023. 

An 8"x8" section of the WFPC2 image is shown in Fig.0 
The DLA galaxy is visible towards the south and the contours 
of the extended Lya nebula detected by PMAS has been over- 
plotted. The exact positioning of the contours was checked 
from the knowledge of the location of the QSO in the PMAS 
images from Fig. ^ The Lya nebula is clearly extended, but 
it is difficult to measure the extension of it in the direction to- 
wards the QSO since it is contaminated by the flux from the 
red wing of the DLA line. Avoiding this contamination requires 
an exact knowledge of the point spread function behavior with 
wavelength such that the flux contribution from the QSO can be 
subtracted. However as argued from Fig. |6l there is no spatial 
overlap between the QSO and the Lya nebula so the contours 
to the north is caused by QSO emission only. In the other direc- 
tion towards the south there is a clear cutoff in the contours. We 
therefore conclude that the Lya nebula is more elongated in the 
east-west direction than in the north-south direction yielding 
the extension of the Lya nebula of 3"x5". Considering that 
the observed extension of the nebula is a convolution of the 
true emission and the seeing during the observations, it could 
originate in smaller clouds in a more complex environment. 

In the WFPC2 image a faint and small galaxy is visible 
at the position indicated in D96 and M0ller et al. (2002), and 
w hich is identical to th e Lyman Break Galaxy (LBG) detected 
by'steiMiliDdilis'). The impact parameter measured in the 
image is 2'.'37 corresponding to 18 kpc in the adopted cosmol- 
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Fig. 6. The left hand panel shows a section of the PMAS spectra consisting of spectra from 100 fibers around the region of the 
DLA. The spectra are oriented with a horizontal dispersion direction. Here are shown 6 groups of spectra where the emission 
from the extended nebula is visible within the small circles. Each group can be considered as representing a Qf.'S wide two- 
dimensional long slit spectrum. In comparison we have in the upper right hand panel shown an artificial long-slit spectrum one 
would have obtained with a I'.'O slit placed in the north-south orientation as indicated in Fig0 The lower right hand panel shows 
the corresponding one-dimensional spectrum from extracting the Lya emission spectrum. This is directly comparable to that in 
D96, in the sense that it exhibits the same red-winged profile, the same line flux, and the same line width. Note that the emission 
within the small circles do not overlap spatially with the QSO emission. 



ogy, the position angle is 159.1°, and the size of the galaxy 
along the major axis is O'.'S, which at z = 3.15 corresponds to 
3.8 kpc. 

In Fig. the size of the Lya nebula appears much 
larger than the size of the DLA galaxy in the WFPC2 im- 
age. To quantify this we investigated whether some low sur- 
face brightness emission is lost in the WFPC2 image due 
to different sensitivities. In the WFPC2 image the 3a lim- 
iting magnitude is 26.8 mag arcsec"^, corresponding to a 
flux limit of 6 X 10~^^ ergcm^^s^^ arcsec^^, while the 
PMAS observations of the Lya object detect emission of 
4xlO~^^ ergcm"^ s~^ arcsec^^, implying that the detection 
limit of the PMAS spectra are roughly the same as for the 
WFPC2 images. Therefore, if the continuum emission were as 
extended as the Lya line emission it would have been detected 
in the HST images. 

The optical magnitude of the LBG corresponds to an L* 
galaxy (D96). In Fig.Qthe DLA galaxy is seen to be composed 
of two components separated by 0'.'25 corresponding to ^2 kpc 
at the redshift of the DLA galaxy, and the two components 
are aligned almost orthogonal to the long axis of the extended 
emission. A galactic outflow will naturally occur along the mi- 



nor axis of the galaxy as described in iHeckman et alJ ( 119901) . 
which would explain the orientation of the elongated Lya neb- 
ula with respect to the LBG orienta tion. The irreg ular morphol- 
ogy of the galaxy was also found bv lMc>ller et al.l (t2002) in their 
HST/STIS data. They noted that the morphology of the DLA 
galaxy was not unusual compared to field galaxies at the same 
redshift. 

Using aperture photometry in IRAF we find the mag- 
nitude of the DLA galaxy mF702W = 24.80 ± 0.1 
using the planetary camera zero point from the WFPC2 
Instrument Handbook. This corresponds to a continuum flux 
of -3x10"^^ ergcm"2 s-i A"^ at A = 1690 A in the rest- 
frame of the DLA galaxy. This is below our detection limit in 
the PMAS spectra and explains the non-detection of underlying 
continuum emission from the DLA galaxy. Dust obscuration 
would make the intrinsic luminosity only larger The derived 
flux from the WFPC2 image, confirming that the UV contin- 
uum is consistent with being flat, is in agreement with the hy- 
pothesis in D96 that the galaxy is very young, and in a star- 
bursting phase. Using the same aperture for the STIS data we 
derive to(V5o) = 25.65 ± O-l in ag reement within Ict of the 
value reported in lMoUer et"al] (l2002l) . 
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Fig. 7. A negative 8"x8" WFPC2 image of the Q2233+131. 
North is up and east is left. The identified absorbing galaxy 
is seen towards the south of the QSO and the contours of the 
object in Fig.^are overlaid. Note again, that the emission is a 
sum of the Lya emission from the nebula and some emission 
from the QSO. The levels of the contours correspond to 5-8(7 
above the median background, with steps of 0.75(7. There is 
a peak in the Lya emission close to the location of the DLA 
galaxy. Indicated by the vertical dashed lines are the position of 
the artificial slit from which the spectrum in Fig.|6lwas created 



6. Lya emission from high redshift objects 

We now compare the properties of other DLA galaxies with 
those for the Q2233+131 DLA galaxy. 

The EW of the observed Lya line is relativ ely large com- 
pared to that of other high redshift galaxies dShaplev et alJ 
I2OO3I) . However, one must note that using a 1" long sUt spec- 
trum one would have obtained a 4 times smaller flux, and there- 
fore decreased the inferred EW by the same amount. Values of 
Lya EWs for high redshift Lya emitters are > 14 A, while 
some may have EWs which are two orders of magnitude larger 
(Kudritz ki et al, 2000). For the Lya emission from the DLA 
galaxies in PKS 0528-250 and Q2206-1958 the measured 
EWs are 63 A and 83 A, respectively, which is typical for LBGs 
at the same redshifts ( M0ner et al 20 02). 

We found that flie hne widfli for the Q2233H-131 DLA 
galaxy suggest larger velocities: 1000 km s~^ compared to 
^700 km s^^ observed for other DLA galaxies. An explana- 
tion for this may be that we have co-added the spectra over a 
large area. Velocity differences of 300 km s^^ have been found 
from one end to the other of the extended emission, and with 
the artificial slit spectra created from the integral field spectra 
we indeed find a smaller value. 

[Leibundgut & Robertson (19991) found evidence for an ex- 
tended Lya emission nebula associated with the DLA in 
Q2059-360 using long-slit spectroscopy, with the slit placed 
at several positions offset from the QSO. They found a velocity 



difference between the Lya emission line and the DLA absorp- 
tion of H-490 km s^^, and that the Lya emission line could be 
described by two components separated by 5 A. However, they 
could not exclude that the Lya emission was affected by the 
QSO, which has a very small redshift difference from the Lya 
emission line. Nevertheless, these observations are remarkably 
similar to those presented here of the Q2233H-131 DLA, except 
for the fact that the QSO in our case is well separated in redshift 
space by Az = 0.15. 

The properties derived for the Q2233H-131 Lya nebula are 
also similar to t he DLA galaxy dete cted in QO 151 -h048 A which 
is described in Fvnbo et al. (7999) and iFvnbo et alJ ( I2000I) . 
They concluded that some of the extended Lya emission could 
be caused by photoionization by the QSO, which has the 
same redshift as the DLA. This is not the case here. The 
proper distance between the QSO and the DLA is 120 Mpc in 
the adopted cosmology. Using the relation between the QSO 
absolute m agnitude and the distance to the DLA given in 
IWarren &'Mc)ller (1996) . we find that the ionizing flux can be 
at most 10~^^ erg cm ~^ s~^ arcsec^-^, which is 4 orders of 
ma gnitude below the o bserve d value. 

IWarren & M0lleil ( 1 19961) found evidence that the DLA 
galaxy termed SI in Q0528-250 has an extension of 1" after 
correcting for atmospheric seeing suggesting that Lya emission 
is more extended than the region emitting continuum radiation. 

Previously a couple of LBGs at 2: — 3.1 ha ve been asso- 
ciated with surrounding Lya emission nebulae in lSteidel et alJ 
(2000), who termed these "Lya blobs". The properties of the 
extended nebula presented here are less extreme, but not very 
different from the Lya blobs. The total luminosity of the Lya 
emission from the DLA nebula is a factor 10 less than from the 
Lya blobs, and the line width of the DLA emission is a fac- 
tor of 2 smaller tohvama et al."2003"). The inferred size of the 
DLA emission is smaller than for the Lya blobs, which have 
sizes of ^ 15", i.e. a factor of 3 larger than the DLA's emis- 
sion, but much of the emission is diffuse, and substructure in 
the blobs is clearly visible. 

It is likely that the observed extension of the DLA galaxy's 
Lya emission becomes larger as one goes to fainter fluxes, and 
since the surface brightness decreases with redshift as (1 + z)^, 
it is difficult to detect faint extended objects especially with a 
4m class telescope. 

7. Origin of thie extended emission 

We address now the question whether the Lya emission is in- 
duced by a superwind from the DLA galaxy and caused by star 
formation, or if it originates in a rotating disk. 

For the assumption of a rotating disk the relation 
Vc = w(FWHM)/(2 sin i) gives = (545 ± 65)/ sin i. 
Such a large value makes the interpretation of a rotating disk 
questionable. A more likely explanation is resonance scatter- 
ing of Lya photons that can also produce a large FWHM of 
the emission line. The resonant nature of the Lya photon nei- 
ther increases nor decreases the total flux emitted by the source, 
only the escape direction, but when the dust-to-gas ratio of the 
environment exceeds 10% of the Galactic value, the extinc- 
tion will be significant due to the increased escape path from 
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the resonance scattering of the Lya photons. While in dust 
free star-forming gala xies the expected Lya EW is 100-200 A 
llCharlot&Falil993h . the fact that the measured EW is similar 
to this value suggests that the most probable source of ioniza- 
tion is massive star formation, and that dust extinction plays a 
small role. 

The red wing of the line profile observed by D96 is 
marginally detected in our lower resolution data. We can repro- 
duce the line profile in the artificial slit spectrum in Fig.|6j but 
it is much weaker when co-adding several spaxels as in Fig.|2l 
indicating the presence of some velocity structure. 

7. /. Lya emission indicative of star formation 

In star-forming galaxies, it is mainly the OB stars that produce 
Lyman continuum photons. These photons will be absorbed by 
the s urrounding neutral hydrogen and re-emitted as line pho- 
tons. |^nnicuta{l99^ gives the relation between the Ha lu- 
minosity (L) of a galaxy and its star-formation rate (SFR). 
Assuming case B recombination (Osterbrock 1989) and that 
the expected ratio of Lya to Ha flux is F(Lya)/F(Ha) « 10, 
implies that 

SFR (Mq yr-^) = L(Lya)/1.26 x 10^^ g^g 

The conversion factor depends strongly on the adopted initial 
mass function, and we adopt an uncertainty of 50%. The rela- 
tion between the SFR and the luminosity also depends on the 
escape fraction of the Lya photons, which is very uncertain as 
the extinction is unknown. 

Assuming that the Lya emission is caused by star forma- 
tion we find SFR ~ 19 ± 10 Mq yr~^, not corrected for 
the unknown escape fraction. This is a typical value com- 
pared to_flie_SFRs found for LBGs having strong Lya emis- 
sion dShaplev et al. 2003.) . It is also consistent with the ob- 
served continuum, as can be shown using the conversion from 
the continuum flux at 1500A in th e restframe of a galaxy 
and the SFR, for which iMadau et all gives the relation 

SFR(Mo yr-i) = 1.3 x lO-^^L^ (erg s'^ Hz"!). We take into 
account an uncertainty for the conversion factor of ^30%. The 
observed R band corresponds to 1500 A in the rest frame of the 
DLA system. Correcting the R band magnitude for Galactic ex- 
tinction and applying an offset converting Vega magnitudes to 
AB magnitudes (F ukugita et al, 1995) cancels each other out, 
i.e. Rab = 24.8 ± 1. Using the relation above this magnitude 
corresponds to a SFR= 12 ± 5 Mq yr~^, which is consistent 
with the SFR derived from the Lya flux. 

Given the uncertainties regarding the conversions from 
fluxes to SFRs, we find that the Lya emission around the DLA 
galaxy could be caused by massive stars. In principle a sig- 
nificant amount of dust would imply that the true Lya flux is 
higher than derived from the UV flux. On the other hand, the 
consistency between the two measurements indicates that dust 
extinction plays no major role in this particular case. A low dust 
content could imply a larger escape fraction of the Lya pho- 
tons. Another explanation of the large Lya lumi nosity could be 
the co mplex dynamics involved in this system jDawson et alJ 
I20n2l) . and not directly related to the dust content as argued for 
nearby galaxies (.Kunth et al..2003>) . 



7.2. Origin of the double peaked emission line 

Several authors have analysed the transfer of L ya photons 
in ne utral clouds (Adams 1972i: IUrhaniak & Wolfe 19*811 
Zhen g & Mira1da-Kscudell2002ir Even in the case of a static 
cloud, the emission profile of Lya would have a double peak, 
whereas an outflow would blend the two components of the 
line. The wavelength difference between the two components 

is AA = \ An, where k is the Boltzmann constant, T 

is the temperature, niH is the Hydrogen mass, c is the speed of 
light, and Aq is the rest frame wavelength of Lya. Assuming a 
temperature of 10^ K, the corresponding splitting of the line is 

0. 2 A at z — 3.15, which is much smaller than the observed 
splitting. On the other hand, given the observed splitting, a 
temperature of 3 x 10^ K is required, which is two orders of 
magnitude higher than temperatures in superwind outflow s at 
distances far away from the central starburst ( Heckman et alJ 
Il990l) . We therefore argue that the observed double peak is not 
caused by a static neutral cloud, but has a more complex origin. 

In the case of an outflow, the double peaked emission line 
can be explained. A galactic wind, powered by star forma- 
tion and supernova explosions, expands outwards and inter- 
acts with the surrounding material, shock heats and acceler- 
ates this, producing emission line filaments. The preferred di- 
rection of expansion is where the pressure gradient is largest, 

1. e. along the rotation axis of the galaxy usually aligned with 
the minor axis. The wind creates a large shell of swept up 
material, and observations towards such an object would in- 
tersect the shell in the front and the rear end giving rise to a 
double peaked emission line profile. Such double peaked pro- 
files of optical emission lines have been observed for nearby 
starburst galaxies which e xhibit these large scal e outflows typ- 
ically termed superwinds jHeckman et aljfl990l) . Additionally, 
a double peaked Lya emission l ine profile was fo und from the 
starburst galaxy T1214-277 bv lMas-Hesse et ab C2003) . who 
concluded that the feature is caused by emission in an out- 
flow. The inferred velocity separation between the two peaks 
from the DLA galaxy emission are similar to the h i ghest ones 
measured for the nearby galaxies in lHeckman et alJ lll990l) and 
similar to the L ya emission line splitting from T1214-277 
jMas-Hesse et al..2003.) . 

7.3. The emitter and the absorber 

In the case of a superwind outflow, the DLA cloud is af- 
fected by its impact, or c ould event be created by a such. 
Taniguchi & Shiova (2001) found that a superwind can create 
Lyman limit systems with Nh > 10^^ cm^^ along filaments. 
Followed by gravitational collapse along such filaments DLA 
clouds can be created. Several authors have pointed out that 
there could be a relation between superwinds from galaxies 
and the metal absorption systems seen in quasar spectra (e.g. 
iHeckman et aD ll990': 'N ulsen et alJll99 8'). thus explaining the 
early chemical enrichmen t of neutral gas cl ouds seen at the 
highest redshifts. However. lBond et all ( 1200 ih conclude the op- 
posite, on the basis of Mg II absorption profiles, that super- 
winds are not causing the majority of DLA systems. 
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The velocity difference between the emitter and the ab- 
sorber could be caused by the impact of a superwind from the 
galaxy which would accelerate a pre-existing low temperature 
clouds in the surrounding environment. Depending on the ini- 
tial distance from the starburst, column density of the neutral 
gas cloud, and the luminosity of the starburst, it can be shown 
that the expe cted acceleration is o f the order of a few times 
100 km s-i llHeckman et alJll990l) . This would be consistent 
with the velocity difference of ^300 km s^^ between the Lya 
emission and the DLA absorption measured here. 

8. Mass estimate 

When the Lya luminosity is known, the mass contained in the 
extended nebula can be eva luated following the method de- 
scribed in lMorse et alJ ( ll998D . The luminosity is given by 

^Lya=jLya^f [ergs"!] (1) 

where JLyQ, is the volume emissivity, V is the volume, and 
/ is the volume filling factor For gas wit h a tempera t ure of 
T — 10^ K assuming case B recombination, IOsterbrockl(ll989l) 

gives 

JLya = 3.56 x IQ-'^^npUe [erg cm-^s"!] (2) 

where Up and are the number densities of protons and elec- 
trons respectively. The total mass of ions in the nebula assum- 
ing negligible metal content is 

Mon = {npirip + nHe-rriHeWf (3) 

Using the approximations rip = Q.9ne and riHe = O.lrip, this 
reduces to 

Mon = l.26nempVf 

« 1.06 X IQ-^'^neVf [Mo] (4) 

Combining the Equations|2and|3]gives 

Afion « 5.91 X 10-46(LLy^V^/)l/2 (5) 

For calculating the volume of the Lya nebula, we make an ap- 
proximation by a cylinder with a radius of 10 kpc and a height 
of 40 kpc. With the calculated luminosity given above we find 
Mion = 5.7 X IQi" Mq. Assuming f ^ IQ-^ taken 

as a ro ugh guess, along with the arguments in lMcCarthv et alJ 
lll990^. we find A/io„ = 1.8 x IG^M©. 

Due to the unknown dust content in the cloud the unex- 
tincted Lya luminosity could be much higher than derived, and 
consequently also the mass. On the other hand, a smaller filling 
factor would decrease the derived mass, i.e. these two unknown 
factors may even compensate each other to some extent. 

9. Conclusions 

Using integral field spectra obtained with PMAS, we have pre- 
sented evidence that the DLA galaxy previously known to be 
responsible for the DLA system in Q2233H-131 has an extended 
envelope of Lya emission. From the spectra we constructed an 
artificial narrow-band image which showed that the extension 



of the emission is 3"x5" corresponding to 23x38 kpc. This 
cloud has a line flux of (2.8 ± 0.3) x 10"^^ erg cm^^ s^^ cor- 
rected for Galactic extinction, corresponding to a luminosity of 
(2.4i°:^) X lO''^ erg s^^ at z = 3.1538. We derive the rest 
frame EW= 190t]l" ^■ 

The source of the ionization is likely star formation within 
the DLA galaxy. The measurement of the Lya luminosity is 
generally considered as a very uncertain method for estimating 
the SFR given the unknown dust obscuration and escape frac- 
tions. With this in mind we find a SFR = 19 ± 10 M© yr~^ 
from the Lya luminosity, a result which is consistent with the 
SFR = 12 ± 5 Mq yr^^ derived from the UV continuum 
flux from the observed ground based R band magnitude. The 
agreement between the two suggests that dust extinction plays 
only a small role. 

A velocity offset of ■~-^270±40 km s~^ between the emis- 
sion and the absorption component in the DLA galaxy does 
not support the hypothesis that the absorbing cloud resides in 
a rotating disk. We investigated the velocity structure of the 
extended Lya emission further by constructing 4 composite 
spectra. By fitting Gaussian profiles to spectra from the eastern 
and western region we found an offset of 2.5 A, which corre- 
sponds to a velocity difference of 150 km s^^. The splitting of 
the spectra into a northern and a southern region gave a veloc- 
ity difference of 300 km s~^. These results combined with the 
extended nature of the Lya nebula is not consistent with the 
properties of a disk similar to that in large spiral galaxies. 

The object's characteristics (luminosity, FWHM, and spa- 
tial extension) are similar to that for some LBGs for which ex- 
tended Lya emission has been detected. These high redshift 
Lya emitters are thought to be created by an outflowing wind. 
Recombination of ionised hydrogen creates Lya photons that 
are resonantly scattered in the environment, which results in 
the observed extended envelope. We therefore consider the in- 
terpretation of an outflow from the galaxy a more likely inter- 
pretation. Furthermore, a galactic outflow can create emission 
lines with double peaked profiles in agreement with the obser- 
vations. 

We argued that the DLA cloud is not a part of a rotat- 
ing disk surrounding the DLA galaxy. The relation between 
the extended Lya emission and the DLA cloud in the line of 
sight towards the QSO is somewhat speculative. The wind re- 
sponsible for the extended emission can create neutral clouds 
with high column densities located along filamentary struc- 
tures. Gravitational collapse in these filaments can create DLA 
clouds. Another explanation could be a previously existing 
neutral gas cloud which has been accelerated by the superwind 
giving rise to a velocity separation of 300 km s^^. A third 
explanation is that the DLA cloud could be an otherwise un- 
related infalling cloud. This third explanation is questionable, 
given the fact that the measured metallicity in the DLA cloud 
suggests some processing of stellar material. 

We have demonstrated the advantages of using integral field 
spectroscopy for investigating DLA systems in the terms of 
confirming the galaxies responsible for the DLA absorption. 
First of all, one can avoid the conventional two step procedure 
for detecting candidate DLA galaxies near the fine of sight 
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towards the QSOs, which later have to be followed up spec- 
troscopically. With integral field spectra one can create artifi- 
cial narrow-band images suited to any wavelengths required. 
Secondly, a non-extended line emitting region could be missed 
by placing the sUt in a less favorable angle, while in the case 
of an extended emission the line fluxes would inevitably be un- 
derestimated. 

We will undertake further investigation of the Q2233+131 
DLA system with PMAS in order to investigate the velocity 
structure of the Lya nebula. Specifically, the system should be 
observed with a higher spectral resolution which should allow 
to determine the presence of P Cygni profiles expected for the 
case of superwind outflows but also in order to investigate the 
double peaked profile in detail. 
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